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ABSTRACT. A probabilistic representation formula for general systems of linear parabolic equa- 
tions, coupled only through the zero-order term, is given. On this basis, an implicit probabilistic 
representation for the vorticity in a 3D viscous fluid (described by the Navier-Stokes equations) 
is carefully analysed, and a theorem of local existence and uniqueness is proved. 



1. INTRODUCTION 

Consider the Navier-Stokes equation in [0, T] x R 3 

V)w + = uAu + f 




(1.1) 

This equation describes, in Eulerian coordinates, the evolution of a viscous incompressible 
Newtonian fluid, where u is the velocity field, p the pressure, / the body force and v > the 
kinematic viscosity. The vorticity field £ = curl u satisfies the equation 

(1.2) d t ( + (u- V)£ = z/A£ + (£- V)u + g 

with g = curl /. As we shall remark later on, the stretching term (£ • V)w can be written in the 
form 

(f ■ V)« = V u ^ 

where V u = |(Vu + Vw T ), which better describes the action of the deformation tensor V u on 
£. The analysis of the vorticity field is a fundamental issue related to questions like the possible 
emergence of singularities (see for instance Beale, Kato and Majda Q, Constantin 0), or the 
description of 3D structures (see for instance Chorin 0). 

The Lagrangian formulation of the fluid dynamics may be important to analyse the vorticity 
field. Strictly speaking, the fluid particles (we mean infinitesimal portions of fluid, not the single 
molecules) move according to the deterministic law 

X(t) = u(t,X(t)). 

However, a virtual Lagrangian dynamic of the particles of the form 

(1.3) dX(t) =u(t,X(t))dt + V2^dW t 

(where W t is an auxiliary 3D Brownian motion) allows us to describe the evolution of quantities 
which are not only transported by the fluid, but have a diffusive character. The vorticity has this 
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property, as many scalars or fields possibly spreading into the fluid. Roughly speaking, we 
prove the representation formula 

at,x) = E[V(t,OMX(0))]+ [ E[V(t,s)g(s,X(s))]ds 

Jo 

where E[-] denotes the mean value with respect to the Wiener measure, £ is the vorticity at 
time zero, X(s) is the solution of equation (11.31) with final condition X(t) = x and V(r, s) is 
the solution of the 3x3 matrix equation 

/ ±V(r,s) = V u (r,X(r))V(r,s), re[s,t] 
\ V(s,s) = I. 

The present paper is devoted to explain the formula in detail, and use it to prove a local-in- 
time existence and uniqueness result. This paper is in a sense the continuation of a paper of one 
of the authors (see Busnello [5|), where the 2D case has been considered. In the 2D case the 
stretching term V u ^ is zero, so V (r, s) = I. The vorticity is purely transported and diffused, 
allowing for a global-in-time control which yields global existence and uniqueness results. In 
Busnello [ 5 1, the probabilistic formula is used to prove such a result, related to the deterministic 
work of Ben-Artzi 01, following a suggestion of M. Friedlin. 

Girsanov transformation is used in a basic part of the work, and the Bismut-Elworthy formula 
is used to treat by probabilistic methods also the Biot-Savart law, which reconstructs u from £ 
(necessary to solve (11.31) ). In the 3D case the Biot-Savart law and its probabilistic representation 
are 

u(t, x) = -±- [ ( g -^ x p) dy= l - T -E[£(*, x + W s ) x W s \ ds. 

In the present paper we extend as much as possible the probabilistic approach of Busnello 
to the 3D case. Now, a priori the stretching mechanism could produce singularities and 
blow-up in £, so we can only work on a time interval [0, r] depending on the size of the data. 
This is the only possible result that also the analytic approaches to equation (11.11) can reach at 
present. Global existence for (11.11) is known only at the level of weak solutions, but we have to 
work at a higher level of regularity to deal with the vorticity. In certain function spaces, global 
existence (and uniqueness) are known for sufficiently small data; in principle the probabilistic 
formulation could lead to such results, but we have found some obstacles, so a probabilistic 
proof of such a result remains an open problem (except for the completely different approach 
of Le Jan and Sznitman lfT9lD . 

The plan of the paper is the following. In Section El we state the precise representation 
formula and the local existence and uniqueness result for the Navier-Stokes equation, with the 
main lines of its proof. However, the full proof of the representation formula and the local result 
are based on three main items that we postpone to the next three sections: 

(i) a general representation formula for linear systems of parabolic equations, given in 
Section |3 

( ii) the probabilistic representation of Biot-Savart law and a number of estimates on it, given 
in SectionlH 

( Hi) a series of estimates for the expected values appearing in the formula for the vorticity, 
given in Section 
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We have chosen this ordering to hi-light the results for the Navier-Stokes equation at the begin- 
ning, for the reader who is not interested in the long list of estimates and preliminaries necessary 
to prove the main theorem. About item ( i) above, we remark that we use a method due to Krylov 
(in the scalar case) that introduces new variables in order to eliminate the zero order terms of 
the parabolic equation. Such method in the case of systems coupled through the zero order 
part is particularly interesting because it reduces the original system to a decoupled one. The 
representation proved in Section|3]can be applied, in principle, to several other systems of equa- 
tions appearing in fluid dynamics, like the equation for u itself (but the term Vp appears in the 
right-hand- side), the equation for the magnetisation variable (see for example [7|), the equation 
for the transport of passive scalars. 

Concerning the literature on the subject, at an advanced state of the work we became aware 
of the interesting papers by Esposito, Marra, Pulvirenti and Sciarretta [10] and by Esposito 
and Pulvirenti ifTTI where a partially similar representation formula was introduced; differently 
with respect to these papers we treat probabilistically also the Biot-Savart law, we use different 
probabilistic tools, we analyse in detail the general case of systems of probabilistic equations 
to understand rigorously the equivalence with the probabilistic representation and we prove the 
local existence and uniqueness result in different function spaces (in particular, for a class of 
less smooth initial conditions). 

There is also a paper by Rapoport [21J dealing with a general class of equations on manifolds 
which in particular throw light on the differential geometric structure of the formula. Also the 
probabilistic representation of systems of parabolic equations has been treated in the literature 
under certain assumptions (our work seems to be more general); see Kahane [ 16 1 and Freidlin 

o. 

Finally, among the literature on probabilistic analysis of PDEs there are possible connections 
with the geometric approach of Gliklikh lfi5l . with recent investigations on vortex method in 
3D by Meleard and co-workers, by Giet (T3, and more closely with a work in preparation 
by Albeverio and Belopolskaya Q where a probabilistic representation for the velocity u is 
employed. Concerning the huge literature on the deterministic analysis of the Navier-Stokes 
equations, results of local existence and uniqueness have been proved in a great amount of 
function spaces, see for instance collections of results in Cannone [6| and von Wahl [ 22 1, or in 
many works of Kato, Solonnikov and many others. We have not found a theorem exactly with 
the spaces used in the present paper, but it may exist somewhere or it may be proved with an 
adaptation of the existing techniques. 

1.1. A physical interpretation of the probabilistic formula for the vorticity. 

1.1.1. Evolution of the vorticity in the non-viscous case. Let us first recall the well-known 
evolution of the vorticity field of an incompressible non-viscous fluid (therefore described by 
the Euler equation). Let x) be the value of the vorticity at time t and point x E R 3 . The 
material point x moves according to the law 

r x(t) = u(t,x(t)) 

1 X(0) = x, 

where u is the velocity field of the fluid. From the Eulerian description of the evolution of £ 
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C(t, X(t)) = V u (t, X(t))£(t, X(t)) + g(t, X(t)) 



(1.5) 



t(t,X(t)) = V(t,0)Mx)+ [ V(t,s)g(s,X(s))ds 



where 



(1.6) 



/ ±V(r, s) = V u (r, X(r))V(r, s) r e [s, t] 

1 V(s,s) = I. 



Take g = for simplicity (the general case is similar); equations (11.41) and d 1 .51) say that the 
initial vorticity £o{%) at point x is transported along the path X(t), and during this motion it is 
modified by the deformation tensor. For instance, the vorticity is stretched when it is sufficiently 
aligned with the expanding directions of V u ; of course the relative position of £ with respect 
to the expanding and contracting (remember that Trace V u = 0) directions of V u changes in 
time, so £(£, X(t)) may undergo a complicate evolution with stretching, rotations, contractions. 
Heuristic reasoning and numerical experiments show a predominance of the stretching mecha- 
nism, and seem to indicate even a blow-up of X(t)) in finite time, for certain initial point 
x. 

If we want to know £(t,x) at a certain time t and point x, we have to solve the backward 
equation 



Jo 

where we have denoted by X t,x {-) the solution of (11.7b . to stress the dependence of the final 
condition x at time t, and by V t,x (r, s) the corresponding solution of equation (11.61) . 

1.1.2. Path integral modification in the viscous case. In the viscous case the position X(t) of 
a material point still evolves under the deterministic equation X(t) = u(t, X(t)). However, the 
vorticity carried by the fluid particle at time t = is not simply transported along its motion 
and modified by the action of the tensor V u ; a diffusion of £ takes place. Let us introduce a 
virtual evolution of fluid particles, subject to u and a random diffusion: 



where W t is 3D Brownian motion. Whether such a motion has a physical meaning or not seems 
to be a similar question to the case of Feynman paths in Feynman integrals approach to quantum 
physics. The initial vorticity ^o(^) decomposes, in a sense, in infinitesimal components along 
the different random solutions of (I1.9I) . proportional to the probability of each evolution (strictly 
speaking such probabilities are zero). If X(t, u) is a path given by (11.91 ). let us denote by P{oS) 




(1.8) 




(1.9) 
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its probability, ignoring for a moment that P(u>) = 0; then an amount of vorticity equal to 
£ (x)P(lo) travels along X(t, u) and is subject to the action V(t, s) of V u along the path: 

at,X(t,u))P(u) = V(t,0)£o(x)P(u)+ I V(t,s)g(s,X(s,u))]dsP(u) 



(the reasoning for the integral effect of g is similar, and we omit it). Now X(t, lo))P(lu) 
is not the total value of the vorticity field at time t and point x = X(t,u), but it is only the 
contribution due to the cj-evolution started from position x: other initial positions and other 
evolutions will reach the point x at time t, and we have to add all these contributions. Therefore 
to compute £(t,x) at a certain time t and point x we have to solve the backward stochastic 
equation 

' dX{t) = u{t,X(t))dt + V2vdW t te[0,t] 
X(t) = x 

to find the various positions X(0,uj) which move to x at time t under different noise paths 
W(t, lo); at the heuristic level each cu gives a contribution x; ui))P(ui) to x) given by 

ri 

rt,x (~L n. , .\ r~ I vt,x fr\. , ,\\ T)l , \ i / t rt,x l~L \ _{ „ \rt,x l 



f (t, x; u))P(u) = V^d 0; ^(^(O; u))P(u) + / V^d s)g(s, X^(s; u)) ds P{lo) 

Jo 

(see (11.81) and (11.41) ). so the total £ (£, x) is given by 

Z(t,x) = E[^(*,0)£ o (^(0))] + [ E[^(t,s)g(s,X^(s))]ds. 

Jo 

This is the heuristic derivation and the physical explanation of the formula. 

2. MAIN RESULT ON THE PROBABILISTIC REPRESENTATION FOR THE VORTICITY 

2.1. Some definitions and notations. First we recall some classical spaces, like the space 
L P (R 3 , R 3 ) of 3D vector fields whose p-power is summable, with norm 

II/IIp=(/ \f(x)\ p dx)K 

</R3 

the space C^'(R 3 , R 3 ) of k-times differentiable vector fields, with norm 

\\g\\ c * = Halloo 
\P\<k 

and finally the space C7^' Q (R 3 , R 3 ) of vector fields whose A; th -order derivatives are Holder- 
continuous with exponent a, with norm 

\\g \\ c k,a = \\g\\ C k + [g]k+a, 

where 

^ \D^(x)-Df 3 g(y)\ 
[9\k+ a = SU P 1, _ ii a • 

|/3|=fc^ eR3 l|X Vl1 

Next we define the spaces where our problem will be set. The velocity field of Navier-Stokes 
equations will be in the space 

(2.1) U a {T) = {ue C([0, T]- Cl(R 3 , R 3 )) n L°°(0, T; C 6 liQ (R 3 , R 3 )) | div«(t) = } , 
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endowed with the norm 



supess < t < T \\u{t)\\ c i, a , 
while the vorticity will be in the space 

(2.2) V a ' p {T) = C([0,T];C fe (R 3 ,R 3 ) nL p (R 3 ,R 3 )) n L°° (0, T; C^(R 3 , R 3 )) , 
endowed with the norm 

sup ess < t < T [|w(t)l|iPncg-, 
where ||'i/ ; ||cf*n£p = \\ip\\p + llV'llc?- We will use also the space 

(2.3) KA T ) = [ueU a {T) supess t < T ||w(t)|| c i,a < M } , 
and the space 

(2.4) VrCO = { e V°' P (T) | supess,< T || ll^ncg- < L } . 

2.2. Probabilistic representation for the vorticity. The formulation of the three dimensional 
Navier-Stokes equations 

d t u - uAu + (tt • V)tt + VP = /, 

(2.5) divu = 0, 

u(0, x) = Uq(x), 
lim u(t, x) = 0, 

|a:|— »oo 

can be given in terms of the vorticity field £ = curl u as 

d t i - uAC + (u ■ V)f - (C • V)« = 
£(0,x) = £o(:e), 
£ = curlii, 
div u = 0, 
lim x) = 0, 

\x\— >oo 

where g = curl /. We shall write the term (£ • V)w as (Vm)^. Moreover, the same term can be 
written as V u ^, where V u is the deformation tensor, the symmetric part of Vu, 

V u = ^{Vu + Vu T ), 

since 

(v«)£ - = ^(v M - v« T )e = e x e = o. 

As we explained intuitively in the introduction (see Section fTTTT) and we shall describe rigorously 
in the sequel, using the representation formula of Theorem I4.4l and the generalised Feynman- 
Kac formula of Theorem 13. 121 the formulation of Navier-Stokes equations can be given in the 
following way: 

(2.6) at, x) = E[a?'*£ (*?■*)] + / nu'fgit - s, X?*)] ds, 

Jo 

i r°° i 

u(t, x) = - / -E[f (t, x + W s ) x W s ] ds, 
* Jo s 
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where the Lagrangian paths (Xf'*) < s <t are processes solutions of the following stochastic 
differential equations 

dXp 1 = -u{t - s, X x /) ds + V2^ dW s , s< t, 

Xq' = x, 

and the deformation matrices (U x,t )o< s < t are the solutions to the following differential equa- 
tions with random coefficients 

dUf = UfV u (t - s, X x /) ds, s < t, 

jjX,t _ j 

u — 1 J 

Here V u is either Vw or the deformation tensor (the name deformation matrices of U x,t refers 
to the latter case). Notice that, with respect to the introduction, we have made a time-reversion 
which simplifies the mathematical analysis. 

A sufficiently regular solution of the classical formulation (12.51) is a solution of (12.61) and 
vice-versa. The main aim of this section is to show that, under suitable conditions, problem 
(12.61) has a unique local in time solution. The claim is proved in the following theorem. 

Theorem 2.1. Given p G [1, § ), a G (0, 1) and T > 0, let £ G Q*(R 3 , R 3 ) n L P (R 3 , R 3 ) and 
g eV a > p {T),andset 

£ o = H&llcPnjy + / lb(s)||cf*ni> ds. 



Then there exists r G (0, T], depending only on Eq, such that there is a unique solution u G 
U a {r), with £ G V a ' p (r), of problem d2~6l 

Proof. We will show that there are suitable L, M and r such that the map 12s o M, where 
& : V2 ,p (r) -> Wft(r) is defined as 

i r°° i 

&(0{t,x) = - / -E[£(t,x + W a )x W s }ds, 
1 Jo s 

and M : W&(t) -> V£' p (r) is defined as 

Afi(tt) (f , x) = E[U!'%(X^)} + f E[Ufg{t - s, )] ds, 

Jo 

is contractive. 

First, in view of Corollary 14.51 M > CL. Using Proposition 15 .51 we see that M maps U^(r) 
toVr(r)if 

(2.7) e 3rM (l + TM)e < L. 

By means of Corollary 14.51 and Proposition l5.61 12s o M is contractive if 

(2.8) CC(is, P )C m (t)s <1, 

where C{y,p) is a constant depending only on p and z/, and lim T ^ Cm(t) = 0. Hence, it is 
sufficient to choose r small enough in order to have both conditions (12.71) and (12.81) verified. □ 

Remark 2.2. As usual, the statement of the above theorem can be read in terms of small initial 
data. More precisely, for each fixed time T, there is a constant e such that if e < e> there exists 
a unique solution u G U a (T), with f G V Q ' P (T), of problem (El) 
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3. The Feynman-Kac formula for a deterministic system of parabolic 

equations 

This section is devoted to the development of a probabilistic representation formula for the 
following system of parabolic equations with final condition: 

(31) f d tv k + \ a ij d l iXj v k + Eti hd Xi v k + {Vv) k + f k = 0, jg = i i 
\ v k (T,x) = ifik(x), x G R d , 

for (i, x) G [0, T] x R d , or the following system of parabolic equations with initial condition 

(3.2) f d t v k = lT,i,j a ij d l i x j v k + J2t=i b id Xi v k + (Vv) k + fk, k = l,. .,1, 
\ v k (0,x) = (p k (x), x G R d , 

where a = aa* and 

a : [0, T]xR d — ► R dxd , b : [0, T] x R d — ► R d , 

(3.3) V : [0, T]xK d — ► K lxl , ip : R d — ► R ; 

/ : [0, r]xR^ R* 

are Borel measurable functions. Additional assumptions will be stated in the sequel. 

At first, for simplicity, assume that / = and all the data are regular. If I — 1, the equation 
(13 .lb - with final condition, has a unique solution given by the Feynman-Kac formula 



v(t,x)=E[<p(xF)e!* v <- r > x r)*-] 
where Xt' x is the solution of the SDE 



(3.4) 



dX** = 6(s, Xl' x ) ds + cr(s, X*<*) dW s , s G [t, T], 



= x, 



where (W t )t>o is a d-dimensional Brownian motion on some filtered probability space. Our aim 
is to extend such formula to the case I > 1. 

Notice that in the case I = 1, for each to, the function 

yt,x _ Q J[ T>{s,X t 3 ' x )ds ^ 

is the solution of the following equation (now V is a scalar) 

1 



(3.5) 



du*'" = ul' x V(r, Xl' x ) dr, r G [t, T] , 

t,x 



So, in the same way, in the case I > 1, we will consider the process [/*»( a, ' y ), solution of the 
equation 

'd[/ r *' (x ' y) = Ur' (x ' Y) V(r, X^)dr r G [*, T], 



(3.6) 



where now both V and £/'>( x ' Y ) are / x / matrices. If Y = I we will write in place of U 1 '^^ . 
Now, the natural conjecture is that, under suitable regularity conditions, the solution of (13.11) is 

(3.7) v(t lX ) = E[U t /<p(X t ^ x )}. 
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In Section im we shall prove (13.71) . under suitable regularity conditions on the coefficients. 
Such formula needs to be modified in order to handle the case / ^ 0, as we show in Section 
13.21 In Section 1331 we shall provide sufficient conditions for the uniqueness of strong solutions 
to system (13 .lb . Finally in Section l3~4l we shall give a Feynman-Kac representation for the 
solutions of the system, with initial condition, (13.21) . 

Remark 3.1. When / = 1, we can write without distinction in formula (13.51) both u l r : x V and 
Du'' 1 , since they are both scalars. If I > 1, the lack of commutativity for the matrix products 
gives that Ul ,x T> and VU*' X are different. The choice in the order of the matrix product in 
equation (I3.6I) . and in formula (13.71) as well, derives from the form of the term V ■ v in system 
<I3 .lb . To have an intuitive idea of this fact, the reader can see the computations in the proof 
of the uniqueness in Proposition 13 .91 (it is convenient to take / = for simplicity). However, 
when one uses backward stochastic equations to represent solutions, the order of U r and T> in 
equation (13.61) changes, see Section l3~4l 

3.1. The homogeneous case. Throughout this section, we will assume 

/ = 

and that the functions b, a and V, given in (13.31) . are Borel measurable functions such that 

(Ai) b, a are sub-linear with respect to x, uniformly in t, 

(A 2 ) b, a are locally Lipschitz-continuous in x, uniformly in t, 

(A 3 ) a is differentiable in x and d Xi a are locally Lipschitz-continuous in x, uniformly in t, 

(At) V is bounded and locally Lipschitz-continuous in x, uniformly in t, 

(A 5 ) (p is bounded and continuous. 

In particular, assumptions (Ai), (A 2 ) and (At) ensure the existence of strong solutions, unique 
in law, for the equations (13.41) and (13.6b . Moreover, from assumption (A 4 ), it easily follows that 

(3.8) II^1r<*< < e T||I,|l ° c , 

where ||£>||oo is the sup-norm. Finally, the previous formula and assumption (A 5 ), imply that 
the function v given by formula (13.7b is well defined and bounded. 
We can now state the main result of this section: 

Theorem 3.2. Assume (Ai)-(A 5 ) and ip E Cb(R d , R 1 ). Then the function 

v(t,x) = E[L/£>(X*f )] 

is continuous and bounded and solves the Kolmogorov equation d3.ll) in the sense of distribu- 
tions, that is 

(3.9) I [ vM*r]dxdt = 0, for all r\ E C£°((0, T) x R d , H l ), 
Jo Jn d 

where 

(3.10) M*rj = -d tV + i Yl €, S"u>l) ~ + 

ij i 

Remark 3.3. The operator M* makes sense since, by assumptions (A 2 ), (A3) and Rademacher 
theorem, the functions 9y ay and <9A are well defined a.e. and essentially bounded in compact 
sets. Moreover, M*r] is bounded in compact sets. 
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To prove Theorem 13.21 we shall use the method of new variables given by Krylov in ifTTl . 
Krylov used such method in order to transform a parabolic equation on R d x [0, T] with poten- 
tial term, into a parabolic equation on H d+2 x [0, T] without potential term. As observed in the 
introduction, we extend this method to systems of parabolic equations. In our case, the elimina- 
tion of the potential term has the additional advantage that the coupling between the equations in 
(13.11) disappears. In other words, we turn system (13.11) into a system of / independent parabolic 
equations on R d +' x/ x [0, T] without the potential term. 

We define the new variables x = (x, Y) £ H d+lxl , and, for each function ip : H d — > R', we 
define the function ^ : R d + 1 * 1 _^ r/ as ^) = Yifj(x). Finally, if u(t, x) : [0, T]xR^ R z , 
we set u(t , x) = Yu(t, x) . 

Prior to the computation of the derivatives of v, we give some notations. We denote by mxn 
the m x n matrix with all entries equal to zero. Given a column vector a £ R d and a I x I 
matrix A, we define the (d + I) (exotic) column vector [%], where the first d rows are given 
by the components of a and the other / rows are the rows of A (the apparent inconsistency 
is inessential, since we shall only use the scalar product defined below). The scalar product 
between two such vectors is defined as 

([%),[%]) = a-l3+(A:B), 
where, as usual, (A : B) = Trace(v4 ■ B) = J2\ ,=i AijBij. 



Given u £ C^^T] x R d ; R') 
du h d{Yu) h 



(3.11) 



9Ya dYij 



since 



d 



^ Y hk u k = 5 ih Uj, h=l,...,l 



it follows that, for each h — 1, . . . , I, the gradient V^ii^ of u h with respect to all its variables is 
given by the following (exotic) column vector 



V x (Yu) h 



V x (Yu) h 

Oix/ 



u 



Oi 



xl 



where the (i-column vector is the gradient with respect to x and the I x I matrix has its rows all 
equal to the /-dimensional vector 0i x ; = (0, . . . , 0) T except for the h th , which is the vector u. 
We want to evaluate next the scalar product ([ ^ v ] , V^Uh) ■ Since 

(YV^dy^Y^h = (YV^jSihUj = (YV) hjUj S lh , 

it follows that 



In particular, if Y 
Let 

(3.12) 



([ Y b v ] ,V„u h ) = b-V x (Yu) h + (YVu) h . 
I, the above quantity is equal to b ■ V x Uh + {Du) 



a(t, x) 



a(t, x) 

0/2 xd 



Odxl 2 

0/2 xt 2 



(3(t,x) 



b(t, x) 
YV(t,x) 



where we understand that a is defined in blocks, where each entry is a matrix itself (notice that 
also D^u h is defined in blocks, and the product (a : D^u h ) is defined as the sum of the four 
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( : )-products of the corresponding blocks). With these notations, if u G C 1,2 ([0, T\ x R d ; R'), 
we have for each h = 1, . . . , I, 

d t Uh + ^ (« : D lu h ) + (P, VxUh) = 

= (Yd tU ) h + a A> ( Yu ^ + E bA > ( Yu ^ + ( YVu ^ 



From this identity it is straightforward to prove that a field u is a strong solution of (13.11) if and 
only if u is a strong solution of system (13.131) . where by strong solution we mean a continuous 
function having continuous first derivatives in time and second derivatives in space, and satis- 
fying the corresponding equation point-wise. In the same way, applying the same ideas used 
above on the adjoint operator, we have the following equivalence. 



Proposition 3.4. A function u is a weak solution of system d3.ll) . with final condition, if and 
only ifu is a weak solution of 

(3.13) d t u h + -(a : D%u h ) + ((3, V w ll h ) = 0, h = l,...,l, 

with final condition u(T,x) = Ytp(x). 

In the sequel we prove that, under suitable conditions, the vector field v(t,x) = Yv(t, x), 
where v is given by (13.71) . is a weak solution of (13.131) . In view of the above lemma, this implies 
that the function given by (13.71) solves system (13.11) in the weak sense. 

The main part is contained in the following proposition, where we relax some regularity 
assumptions on the coefficients of a theorem of Krylov iHTl . Indeed, the drift and the diffusion 
defined in formulae (13.121) are neither bounded nor globally Lipschitz-continuous, in contrast to 
the assumptions of ifTTI . The same problem occurs for the final condition. On the other hand, 
both the drift and the diffusion are locally Lipschitz-continuous and with linear growth (in all 
variables, including Y). 

Proposition 3.5. Let m G N and consider the scalar parabolic equation 

(3.14) d t u + -(a : D 2 u) + (J3, Vu) = (t, x) G [0, T] x R m 

with final condition u(T, x) = ip{x), where a = 77* and 

f3 : [0, T] x R m -> R m , 7 : [0, T] x R" -> R mxm , ^ : [0, T] x R m — > R, 
and assume that 

(i) (3, 7 are Borel measurable, sub-linear and locally Lipschitz continuous in x, uniformly 
in t, 

(ii) ip is continuous and with polynomial growth, 

(iii) j(t, •) is continuously differentiate for each t and d Xi ^ are locally Lipschitz continuous 
in x, uniformly in t. 
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Set u(t, x) = ~E[ip(Z^, x )], where Z l / X is the solution of the SDE 

' dZ l r > x = 0(r, Zl: x ) dr + 7 (r, Z l ; x ) dW r , re [t, T] , 

ryt,X 

Z/ f — X, 

where (Wt)t>o is an m- dimensional standard Brownian motion. Then u is a weak solution of 
(EH; for each r, e C c °° ((0, T) x R m ), we have 

r-T p 

uN*r] dx dt = 0, 



JH 



N*rj = -d tV + \ dl Xj (a ljV ) - dxtiM- 





X 


1° 


\x\ 



where 

Y ' // = —t), ii - ^ 

Proof. If everywhere in the assumptions of the proposition we have global Lipschitz-continuity 
(instead of local Lipschitz-continuity), the proposition follows from Theorem 5.13 of Krylov 
ilTTl . In the general case, we proceed by truncation. Let \& n e C°°(R m ) be such that 

and set (5^ = ty n f3 and 7^ n ^ = ^n). Fix a Brownian motion (fi, T , Tt, W t , P) and denote 
by Zl' x ' n the solutions to the corresponding SDEs. The sequence Zp x ' n converges to Z^' x in 
probability uniformly on compact subsets of [0, T] x R m . 

Suppose first that ifi is bounded. Then u n (t,x) = E,[ijj(Zp x ' n )} converges to u(t,x) and 
Pip converges to f3 Xi , d Xi a^ to d Xi a and d XuXj a^ to d x . jXj a uniformly on compact subsets of 
[0, T] x R m . Let rj e C^°, since N*r) is abounded sequence (see Remark l33l) . by the dominated 
convergence theorem, J u n N*r) converges to J uN*rj, where N* is the operator corresponding 
to the approximate coefficients. Since u n are weak solutions, it follows that u is also a weak 
solution. 

If ip is not bounded, we take a sequence of bounded continuous functions ijj n — » ip such that 
|^n(^)| < IV'C^)!- From Theorem 4.6 of Krylov ifTTl . we have E[|Z^ x '| fc ] < c(l + \x\ k ), so that 
u n {t, x) < c(l + \x\ k ) by assumption (u), and again we conclude by the dominated convergence 
theorem. □ 

We are now ready to prove the main theorem. 

Proof of Theorem \3?2\ First we show that v is bounded and continuous. The boundedness 
comes from (13.81) and the assumptions on cp. In order to show the continuity, we take a se- 
quence (x n , t n ) converging to (x, t). From Lemma 2.9 of Krylov ifTvl . the function (i, x) — ► 

(X l ' x , U t,(x,I) ) e C([0,T},R d+lxl ) (where by convention (X^, ui Ax,I) ) = (x,I) if a < t) is 
continuous in probability. Hence, there is a subsequence such that convergence is almost sure. 
Finally, the conclusion follows from the bound (13.81) . the assumptions on ip and the dominated 
convergence theorem. 

We show then that v is a weak solution. We have the following two ingredients: 

(i) the two systems of SDEs (13.41) and (13.61) can be thought as a unique system where 
the solution (X*' x , Ur takes values in Jl d+lxl and drift and diffusion are given by 
formulae (l3~T2b . 
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(ii) Since by uniqueness for equation (13.61) it follows that U^'^ = YU^ X , for the function 
v defined in (13.71) . we have 

v{t,x) = Yv{t,x) = E[YU^<p(Xt X )) = E[U^ B ' Y) <p(X^ B )) = E^X^, ufi x ' Y) )]. 

From these two facts, by Proposition l3.51 v is a weak solution to system (13 .13b . By Proposition 
3.41 v is a weak solution to system (13.11) . □ 

The regularity assumption (A A ) on the term V can be relaxed with the following condition 
(A' 4 ) V is bounded and uniformly continuous. 
In fact we can deduce the following corollary. 

Corollary 3.6. Assume (Ai)-(A 3 ), {A' A ) and (A$). Then the function 

v(t,x) = EfE^VGXr 5 )] 
is continuous and bounded and solves the Kolmogorov equation d3.ll) in the sense of distribu- 
tions: 

[ [ vM*r]dxdt = for all rj G C~((0,T) x R d ,R d ). 
Jo Jn d 

Proof. Let p n be a sequence of mollifiers and set V n = V * p n and v n (t, x) = E[U^ x n if(X^, x )], 
where is the solution of (13.61) corresponding to V n . 

Since V n — > V uniformly in [0, T] x H d , we have U t '* — ► C/ t ,a: in L 1 (fi), uniformly in 
[0, T] x R d . Consequently, v n (t, x) — > i;(£, x) and P n u n — >■ uniformly [0, T] x H d . Since v n 
are weak solutions of the corresponding approximate problem, in the limit v is a weak solution 
of Mv = 0. □ 



3.2. The inhomogeneous case. In this section, Theorem l3.2l wiil be extended to the inhomoge- 
neous case. We will show a Feynman-Kac representation formula for the complete system (I3.lt . 
that is with / ^ 0, with final condition. Throughout this section we will assume (Ai)-(A 3 ), 
(A' 4 ), (A 5 ) and the following 

(A e ) f : [0, T] x R d — > H l is bounded and uniformly continuous. 
Theorem 3.7. Assume (Ai)-(A 3 ), (A' 4 ), (A 5 )-(A 6 ). Then the function 

(3.15) v(t, x) = E[U^<f(X^)} + ^ Ep^fir, X?)] dr 
is a weak solution of (13.21) . that is, 

[ [ (uM*r] + frfjdtdx = 0, r] G C£°((0,T) x R d ,R ; ). 

The main idea to prove the theorem is to introduce a new component (we apply again the 
method of new variables of Krylov ifTTlO and prove that v is a solution of system (13.11) if and 
only if v — (v x , . . . , 1) solves the system 

(3.16) d t -t; + - ^ <!,,(),. , r /r + ^ hd Xi v + (X>v) fc = 0, 



2 



with final condition v(T, •) = (</?i, . . . , fi, 1), where Z> = ( g J). Notice that = ( ^ +/ ) , so 
that the component vi + \ is obviously a solution. 
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The key lemma is the following. 

Lemma 3.8. The function v = (vi, i>2, vi, 1) is a weak solution of (13.161) if and only if 
v = (vi, i> 2 , .:,Vi) is a weak solution of (I3.ll) . 

Proof. A weak solution of (13.11) is a function v such that J J (vM*rj + frj) = for each test 
function 77, or equivalently J f (vL*rj + vV>*r\ + frj) = 0, where the operator M* has been 
defined in (13.101) and L* is defined as 

l* v = -d t r] + - d liXi K 7 ?) - 9 *i( 6 ^)- 

Let fj = (rj, rji + i) be a R z+1 -valued test function. Since V* = ( ^* ° ) , we have 

vV rj = I ' Q 

It comes out that t> is a solution of the inhomogeneous equation if and only if v solves J f (vL*rj+ 
ffj) = 0, that is, if and only if v is a weak solution of system (I3.16I) . □ 

We can now prove the main theorem of this section. 

Proof of Theorem 1X71 Let (p be the function (cpi, ipi, 1) and U l s ' x be the solution of 



(3.17) 



d fjt,x = [/t,*p( Sj x t,^ dS) s G ^ T ] ; 



Since X? and / satisfy assumptions (A 4 ), (A 5 ) and (A 6 ), the functions and X> satisfy 
assumptions (A 4 ) and (A5). Hence, by Corollary 13 .61 the function 



(x,t)-E[^(Aj B )] 



is a weak solution of system (13.161) . 
We write Ul' x in blocks: 

(/It, IE 
°s 

where A s is a Z x I matrix, b s G R d is a column vector, c s G R d is a row vector and d s is a 
scalar. With this position, the Cauchy problem (13.171) is equivalent to 

( dMf = A*fV(8, Xl' x ) ds, Af = I h 

Ut,x _ At,xff n vt,x\ J„ \t, x 



and it is easy to see that 



db t f = A t > x f(s,X t s ' x )ds, &r = 0, 

ddf = cl' x V(s, Xl' x ) ds, cf = 0, 

ddF = <*ff(s,X¥ ! )ds, df = \, 



A t,x = jjt,x b t,x = J upoffa x y x ) dr 



n 

c tx = fjtf = L 
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Consequently, 

E[U^<f(X^)} = E[( ut f 6 f )(^ x p)] =E[( u T X ^ xt f+ b T X ))] 

□ 

3.3. A uniqueness result. In the preceding sections, we were concerned with the existence of 
a weak solution of the parabolic system (13.11) having a nice probabilistic representation. The 
aim of the present section is to provide sufficient conditions for the uniqueness of solutions. 
In Proposition 13.91 we shall see that the strong solution, if exists, is given by our probabilistic 
representation, hence is unique. In Theorem 13 .101 we will show, under some special conditions 
on the coefficients, that weak solutions are also unique and are given by the probabilistic rep- 
resentation. Such special conditions on the coefficients are satisfied in the application of the 
probabilistic representation to the Navier-Stokes system: if the velocity field is regular enough, 
the coefficients in the equations for the vorticity satisfy the special conditions. Hence, for each 
fixed regular velocity, there exists a unique weak solution of the vorticity equation given by the 
Feynman-Kac formula. 

Let C 6 1,2 ([0, T] x R d , R) be the space of continuous functions having first and second deriva- 
tives in x and first derivative in t continuous and bounded. We start by showing that, if the 
solution of the parabolic system is regular, then it is given by formula (13.151) . 

Proposition 3.9. Let v e C b 1,2 ([0, T] x R d , R) be a strong solution of system (13.11) . with final 
condition. Then v is given by formula (13.151) . 

Proof. It is sufficient to show that the process 

U*?v(r,X<?) + J' 'u?f(8,X?)ds, re[t,T], 

is a martingale. Indeed, if h G {1, ...,/}, by Ito formula (we omit for simplicity (r, X^: x ) from 
the term v(r, X*' x ) and from the coefficients, and the subscript r from the term Ul' x ) 

d r (U^v) h = Y, d ^ U ^) = Y,^dvk + v k dU^) 

k k 

k i i,j 

dWi]+Y,v k Ut x V ik dr 

i,j k,i 

= ~ J2 U hk (/* + VkiVl ) dr + ( dM ^ h + Yl v k U hi V ik dr 

k i k,i 

= -d r (J\lJ t > x f) h ds) + {dM r ) h 

since v is a solution of system (13.11) : (M r ) re u,T] is me (i-dimensional martingale, vanishing at 
r — t, given by 

(dM r ) h = £ d Xi v k a l3 dWl 

k i,j 
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Moreover, M r is square-integrable, since v G C^' 2 , Uj. x is bounded by (13.81) and 



t,x\2] 



sup E[|X, 

t<r<T 

is bounded. □ 

Theorem 3.10. Let ip be bounded and continuous, f and T> bounded and uniformly continuous. 
Suppose that a is constant and b is a Borel measurable and Lipschitz-continuous in x function 
such that div b = 0. Then the function 

v(t, x) = E[V**(p{XF)] + f E[E£*/(r, X?)] dr 
is the unique weak solution of the parabolic system (13 . 1 b . 



The proof of the theorem is based on a regularisation by convolution, in order to apply the 
uniqueness result of the previous proposition. 

Let p G C°°(R d , R), < p < 1, with support in the ball of radius one, and set p n (x) = 
n d p(nx). Let J n be the convolution operator: J n (u) = p n * u. 

Lemma 3.11. Let b : R d — > H, d be a Lipschitz-continuous function, such that div b = (in the 
sense of distributions). Then there is a constant C such that for each u G Cb(R d , R/), 

(3.18) \([J n ,b- V]u)(x)| < C sup \u{y)\ foralln, 

y eB 1/n (x) 

where [J n , b ■ V]m = J n ((b ■ V)m) — {b • V)J n u is the commutator. Moreover 

(3.19) [J n , b • V]w n — ^ uniformly on compact sets 
Proof. Fix u G Cb(R d , R'). Since div 6 = 0, by integration by parts we have 

([J n ,b-V]u))(x) = 

= {pn *{b-V)u-{b-V){ Pn *u))) (X) 

Pn(x - y)(b(y) ■ V y )u(y) - (b(x) ■ V x )(p n (x - y))u(y) dy 
V y pn(x ~ y)(b(x) - b(y))u(y) dy. 

R d 

Taking the norms in R' we get 

\([J n ,b-V]u)(x)\ < [ \V Pn (x-y)\-\b(y)-b(x)\-\u(y)\dy 

JB 1/n (x) 

< cLUVplloo sup \u(y)\ 

yeB 1/n (x) 

where L is the Lipschitz constant of b. So far, we have proved (13.181) . Concerning (13 . 1 9b . it 
is easy to see that the claim is true for u G C£°(H d , H l ). If u is only C&, the claim follows 
from approximation with functions (in the sup-norm, on compact sets) and from the bound 
(l3~T8l) . □ 

We apply now the previous lemma to prove the main theorem. 
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Proof of Theorem \3. 1 0\ Let v be a bounded and continuous weak solution of system (13 .lb . The 
sequence v n = p n * v belongs to C([0, T], C£°(R d , R)) and v n — > v uniformly on compact 
sets. We want to show that v n is a weak solution of 

(3.20) d t V n + - ^ a ij9l iXj V n + ^ b $*i V n + VVn + Pn* f + Rn = 0, 

with final condition v n (T) = p n * ip, where R n = [J n , b ■ V]i> + [J n , T>)v. Indeed, v is a weak 
solution of d3.lt . so that we can use Q n = p n * r\ as a test function, where 77 is again a test 
function and p n (x) = p n (—x), to obtain with some easy computations 

r f (vM*c+fo= 

JO JR, d 

*ii » 

i,3 i 
+ ff V(jnf+ [Jn, b-V]v+ [J n , V)v) 



(notice that for each u, J u(p n * 77) — J rj(p n * u)). 

Since v n belongs to C([0, T], C£°(R d , H 1 )) and p n * f + R n is bounded and continuous, we 
argue that the distributional derivative d t v n is bounded and continuous and, therefore, a strong 
derivative. Hence v n G C\ ' 2 and it is a strong solution of (13 .20I) . Proposition l3.9l vields 

V n (t,x) = E[U$° Pn * ip{X*f)\ + y nUr X {Pn * f + Rn)(Xf?)] dr. 

It is easy to check that \\[J n , D]v\\oo < 2||£J|| 00 ||t;|| 00 and [J n ,D]v — > 0, uniformly on compact 
sets. Hence, by the previous lemma, i? n is bounded, independently of n, and i? n — > uniformly 
on compact sets. Using (I3.8I> and the dominated convergence theorem, we obtain 

v(t,x)= lim Vn (t,x) = E[4'V(^)]+ f nuTf{r,Xl' x )]dr. 

□ 

3 .4. The formula for parabolic systems with initial condition. In this section we describe the 
probabilistic representation of weak solutions to the system (13.2b . with initial condition. Indeed, 
in the sequel we will use the results of this sections to give a probabilistic representation for the 
solutions to the Navier-Stokes equations, which is a parabolic equation with initial condition. 

We will obtain the representation formula for the forward parabolic system using the repre- 
sentation for the backward parabolic system and a time inversion of the coefficients. To this 
aim, we will consider the following stochastic differential equations 



(3.21) 



dXp^ = b{t - r, X^) dr + a[t - r, Xp**) dW r , re [s, t] , 
x ,x,t = X) 
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and 

rS,(x,Y),t t rS,(x,Y),ts 



(3.22) 



jjS,{x,Y),t y 



where, as usual, JJ 8 ^^' 1 = U s,x,t when Y = I. 

Theorem 3.12. Let the data b, o, ip, V and f satisfy assumptions (Ai)-(A 3 ), (A' 4 ) (in paseUH 
(A 5 ) (in page^ and (A 6 )(in paseUH). Then the function 

(3.23) v(t, x) = E^'VC*?**)] + f E[U^f(t - r, X^)) dr 

Jo 

is a weak solution of (13.21) . with initial condition. 

Moreover, if a is constant and b is globally Lipschitz-continuous in x, then v is the unique 
weak solution. 

Proof. Let v(t, x) = v(T — t, x). If v is a weak solution of (13.21) . by easy computations it follows 
that v is a weak solution of 

d t v(t, x) + \ V ciijiT - t, x)d 2 x v(t, x) + V bi(T - t, x)d Xi v 

(3.24) i 

+ V(T - t, x)v(t, x) + f(T - t, x) = 0, 

for t e [0, T], with final condition v(T, x) = <p(x) (and vice-versa). 

By Theorem l3.71 a solution v of (13.241) is given by the following formula 

v(t, x) = E[U¥' T ip(X¥> T )} + ^ E[U^ T f(T - r, X^ T )} dr, 

where Ul ,x,T and X^ X ' T are given respectively in (13.221) and (I3.21I) . We can conclude that a 
solution v of the forward parabolic equation (13.21) is given by 

v(t, x) = E[U^ T ^(X^ T )} + f E[Uj-^ T f(T - r, Xj~^ T )] dr 

JT-t 

Finally, one can easily check that, for each r G [T — t,T], the joint law of the random variables 
rjT-t,x,T an( j xJ~^ x ' T is equal to the joint law of the random variables U®+f_ T and X®+f_ T . In 
conclusion, formula (13.231) holds. □ 

The representation formula above appears more complicated than the formula for parabolic 
systems with final condition (13.151) : the stochastic processes X r in (13.151) are the solutions of 
a fixed SDE corresponding to different initial conditions, while the stochastic processes X r c ' x,t 
and U®: x,t in (13.231) solve for each t a different SDE. A different representation can be given, 
which is more appealing at the heuristic level, even if less suitable for stochastic calculus. 

Consider the following backward SDE 

(3.25) Yf' x = x + f b(s, Y*' x ) ds + f a(s, Y* ,x ) dW s , re[0,t], 

J r J r 

where dW s denotes the backward stochastic integral with respect to the Brownian motion W s 
(see Kunita [ 18) for the definition of the backward integral). Notice that the final condition 
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Y t ' = x has been imposed here. Let V^" x , 0<s<r<tbe the solution of 
(3.26) 



d ys,t,x = V ^ yt^yate rf r, r G [s, t] , 

ys,t,x j 



Theorem 3.13. Under the same assumptions of the previous theorem, a weak solution of the 
parabolic system (13 ,2b . with initial condition, is given by the following formula 

(3.27) v(t, x) = EiV^'MY^)} + f m r>t ' X f(r, Y^)] dr, 

Jo 

where Y t)X and V r,t ' x are given respectively by (13.251) and (13.261) . 

Remark 3.14. We want to give an interpretation of the representation formula given above. 
Suppose for clarity that / = 0. Consider the trajectory Y^' x {uj) of a virtual particle which is in x 
at time t, transported by a velocity field and subject to a diffusion, and evaluate v(Q, Yq X {uj)) = 
(P(Yq X (uj)). Then we take into account, through the vector field V^' t,x , the effects of the tensor 
V along the given trajectory in the time interval [0, t). Finally, by taking the expectation, we 
consider the mean effect of all virtual particles. 

Before giving the proof of the theorem, we need the following simple lemma for the time 
inversion of a stochastic integral. 

Lemma 3.15. Let (W s ) s >o be a Brownian motion. Fix t > and set 

B s = w t -w t ^ s se[o,t\. 

Let TY = cr (W r | r G [0, s]) and Tf t = o [B u — B v \ s < v < u < t) and let g(s) be a con- 
tinuous and bounded process adapted to the filtration F^. Then the process f(s) = g(t — s), 
s G [0, t] is Tf t -adapted and for all a, b such that < a < b < t, 

g{s)dW s = [ f(s)dB s . 



Jt-b 

Proof. Since B u — B v = W t - V — Wt~ u , we have T^_ s = Tf t and this gives the first statement. 
Take now a sequence of partitions of the interval [a, b\. 



TV r 



:{a = s n <s n 1 <...<sl=b} 



such that |7rJ — > 0. We have 



[ g(s)dW s = \imJ29(s?)(W s?+1 -W s n) 

J a 

= lim Y,9(t ~ r?)(Wt-r f+1 ~ W t . r? ) 

n— >oo ' * 1+1 1 

lim J2M)(B r? -B r ^) 
f(s)dB„ 



n—foo ■ 
-t—a 



t-b 

where rf = t — sf , i — 1 . . . k n . □ 
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Proof of Theorem \3.13\ We need only to show that 

X°%* = Y?* and U^f = V[' X '\ P - a.s. 

since such formulas, formula (13.231) and a change of variables give us (13.27b . 

We prove the first equality. Fix a Brownian motion (W r ) r >o and consider the solution X®: x,t 
of equation (13 .21b . By Lemma (13.151) above, it follows that X^f/ satisfies the backward SDE 
(13.251) with respect to the Brownian motion B s defined in Lemma 13.151 Since equation (13.251) 
has a unique strong solution, we have the first equality. 

We proceed to prove the second equality. Fix u so that r — > Y^' x (uj) is continuous. The key 
observation is that 

V r s ^ x {to) = V^ x {lo){V^ x (uj))-\ 0<s<r<t, 

and it is true since 

with initial condition (V° ,t,x (u)y 1 = I, so that it easy to check that V°' t ' x (u)(V^ t ' x (u))-' 1 
satisfy equation (13.261) . Finally, by evaluating 

d P K^(w) = dr[^(w)W^(w))- 1 ], 

we see that both V^ x {u) and r -> U^;\uj) solves the ODE: 

dU r = -U r V(r, Y*> x (uj)) dr re [0, t], 
with final condition U* = I. □ 



4. a probabilistic representation for the newtonian potential and the 

Biot-Savart law 

In the present section we aim to give a probabilistic representation for the velocity field of an 
incompressible fluid in terms of the vorticity field £ = curl u. 

Under suitable assumptions on £, the Poisson equation —Atp = £ has a solution, given by 



4tt J r3 |x - y| 

(ip is a vector field and the equation is interpreted component- wise). Let u(x) be defined as 

u(x) = cmlip(x), i.e. 

If div£ = 0, then also div^> = and divu = 0, and this implies also curl curl?/) = —Aip. 
Therefore curl u = £, i.e. u is the divergence-free velocity field associated to £. The equality 
(|4~TT) is the Biot-Savart law. 

In order to give a probabilistic representation of this formula, it is necessary to give a repre- 
sentation of the solution of the Poisson equation and of its derivatives. 
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4.1. A probabilistic representation for the Newtonian potential. In this section we study a 
probabilistic representation of the solution of the Poisson equation. The deterministic regularity 
results are classical (see for example Gilbarg and Trudinger lfT4ll and Ziemer ll23l0 . so we will 
focus on the probabilistic formula. 

Let / : R 3 — > R be an integrable function. We define the Newtonian potential with density 
/as 



Nf(x) = j- [ _J_/(y)dy. 

4tt J r3 \x - y\ 



If / is regular and with compact support, Nf is a solution of the Poisson equation. 

Let A = |A, it is well known that A generates, on the space C (R 3 ) of all continuous 
functions vanishing at infinity, the strongly continuous semigroup 

P t f(x) = E[f(x + W t )] x G R 3 , t > 0, / G 0,(R 3 ), 

where (W t )t>o is a 3D-standard Brownian motion. The resolvent of A can be written as 

POO 

((A-\I)- 1 f)(x)= / e- Xt E[f(x + W t )]dt, /gC (R 3 ), 
Jo 

so that we can argue that the integral 

POO 

(4.2) / E[f(x + W t )]dt. 



converges to A~ l f(x) = 2Nf(x) (indeed, at this stage, we do not know if A is invertible). 

As a first step, we find some conditions on / in such a way that formula (14.21) produces a 
solution of the Poisson equation. 

Proposition 4.1. Let f e L P (R 3 ) n L q (R 3 ), with l<p<|<g<oo. Then the integral in 
(14.21) 15 convergent for all x G R 3 and is equal to 2Nf(x). Moreover Nf G Cq(R 3 ) and 



\\Nf\\ 00 <C m (\\f\\ p + 
Proof. For every r > 1, by Holder inequality, 

(4.3) E\f(x + W t )\ = — -i™ / |/(x + y)|e-^l 2 ^<C7 r r^| 

so that, by using the above inequality with r = p and r = q and by integrating by time, 

POO pi POO 

E\f(x + W t )\dt< E\f(x + W t )\dt+ E\f(x + W t )\dt<C(\\f\\ p +\\f\\ q ). 
Jo Jo Jl 

This will prove also the final inequality, once the other properties are verified. The integral in 
(14.21) is equal to 2Nf(x) since (we can use Fubini theorem because of the previous inequality) 

/OO P POO i 

E[f(x + w t )]dt] = J n3 f( x +y)J o j^t]^ e ~ ltlyl2dtd y 

J—f(x + y) dy = 2Nf(x). 
R3 2n\y\ 
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We know from Gilbarg and Trudinger lfl4ll . that / G L g (H 3 ) implies, by Sobolev embeddings, 
Nf G C(R 3 ). The behaviour at infinity is less standard, so we give a probabilistic proof of it. 
Thus let us show that Nf G C (R 3 ). Indeed, for each R > 0, 

poo poo poo 

I E[f(x + W t )]dt= / Ef{x + W t )I {m>R} dt+ / Ef{x + W t )I {m < R} dt 
o Jo Jo 

and, in order to show that Nf(x) converges to as |se| — > oo, we will prove that the first term 

converges to 0, uniformly in x, as R — > oo, and the second term converges to as \x\ — ► oo for 

each R > 0. 

For the first term the claim is true since, as in (I4.3I) . 

sup E\f(x + W t )\I {lWtl>R} < C{\\f\\ P + \\f\\ q )(t- 3/2p I[x,oo)(t) +r 3/a9 / [0 ,i)(«)) 
xeR 3 

and 

sup E\f(x + W t )\I {lWt]>R} < Cr 3 / 2 ||/|| p ( / e-kW) W ^0 
xeR 3 J\y\>R 
as R — > oo. As regards the second term, we can proceed as in (14.31) and bound the term 

E|/(x + W t )\I { \ Wt \<R} with 

C ^ 1 1/(2/) ^-[ly-^l <«} I Ip^i.oo) (*) +t~^\\f(y)I{\ v - x \<R}\\ q I[o,i)(t)), 

so that, after the integration in time, the above term converges to 0, since / G L P (R 3 ) n L q (R 3 ). 

□ 

In the second step, we study the derivatives of Nf. Notice that, for a regular /, Bismut- 
Elworthy formula (see for example 0) gives 

D Xi E[f(x + W t )} = ~E[f{x + W t )Wl\. 

In this simple case, with the Brownian motion, such formula can be easily checked by means of 
the Gaussian density. 

As in the previous proposition, one could expect that, under suitable conditions, it is possible 
to write the derivatives of Nf with the probabilistic representation suggested by the formula 
above. Indeed, this is the case, as the following proposition shows. 

Proposition 4.2. Let f G L P (R 3 ) n L q (R 3 ) for some l<p<|<3<g< +oo. Then 
VNf G C (R 3 ) and for each x G R 3 , 

(4.4) 2D Xi Nf(x) = / -E[f(x + W t )WH i = 1, 2, 3. 

Jo 1 

Moreover 

(4-5) l|ViV/||oo<C M (||/|| p +||/|| ff ) 
Proof. By Holder inequality, 



jE\f(x + W t )W*\ = ^ I f(x + y)y l e 



£M 2 



R3 



(4-6) < ^\\f\\ p Vit% 

< C\\f\\ p t-i-* 
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and, as in the proof of the previous proposition, the time integral is finite and bounded with 
respect to x, by the assumptions on p and q. Moreover it can be easily seen, by the same 
arguments used in the previous proposition, that formula (14.41) and inequality (14.51) hold and 
that ViV/ G C (R 3 ). " □ 

In the last step, we study the second derivatives of the Newtonian potential. The regularity of 
the following theorem is based on the classical Schauder estimates. 

Proposition 4.3. Let f G L P (R 3 ) n Q*(R 3 ), with 1 < p < §. Then Nf G Cf '"(R 3 ) n C (R 3 ), 

||iV/|lc^(R3)<5(||/|Up(R3) + ||/||c r (R3)) 

and Nf is the unique solution of the Poisson equation in Co(R 3 ) H C7 2 (R 3 ). 

Proof. From the previous proposition, we know that Nf G C^(R 3 ). Bismut-Elworthy formula 
gives us 

(4.7) D XiX .Ef(x + W t ) = ^E[(D x ^)(x + Wt)W{], 

t 2 2 

where ip(x) = E/(x + Wt). Hence, in order to show that 

f°° 1 

D XiXj Nf(x) = / -E[(D x .iP)(x + Wt_)Wl\ dt 
Jo f 22 

holds, it is sufficient to show that (14.71) is integrable in time in the interval [0, oo). 

First, by the Bismut-Elworthy formula, we see that 

D ai iP(x) = *E[f(x + W*)Wi] 

and, by (|4.6I> . that 

(4.8) H^^lloo < Ct-^-^H/Hjp. 
Moreover, since / G C7^(R 3 ), 

\D x ^(y) - D Xi i;(x)\ = ~E\f(y+W,)-f{x + Wt)\-\Wi\ 

(4.9) < Cr^[f] a \x-y\ a . 
Now we show that (14.71) is integrable in time. By (14.81) 

~E\(D Xi ip)(x + Wt)wi\ ^cr^WfWpEiwii^cr^WfWp 

t 2 2 2 

and (I4.7I) is integrable in [1, oo). By (14.91 ) it follows that 

^E\(D Xi i;)(x+Wt)W{\ = ^E\[(D Xi iP)(x + W*)-(D Xi iP)(x)]Wi\ 



t 



< criE\Wt\ a \wi 



< cr 1+ ^[f] a 

and (14.71) is integrable in [0, 1). 

In conclusion, the probabilistic representation formula for the second derivatives holds and 

WD^Nfw^Kcdifw.+ifu). 
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By Schauder's theory, since Nf G Cf(R 3 ) and / G C b a (R 3 ), it follows that Nf G Cf a (R 3 ) 
and 

ll^/llc^(R3)<C(ll/ll P +ll/llc f ) 

(see for example Lunardi lEUl ). Moreover, Nf solves the Poisson equation (Lemma 4.2 of 
Gilbarg and Trudinger ifbHO and the solution is unique by the maximum principle. □ 

4.2. A probabilistic representation for the Biot-Savart law. We apply now the theory devel- 
oped in the previous section. The following theorem, which is actually a mere corollary of the 
above results, is nothing but the well known Biot-Savart law. 

Theorem 4.4. Let £ G L P (R 3 , R 3 ) n Q*(R 3 , R 3 ), with 1 < p < § and < a < 1. There is a 
unique u G C b 1,a (R 3 , R 3 ) n C (R 3 , R 3 ) such that 

curl u = £, div u = 

and such solution is given by the following formula 

1 f°° 1 

u(x) = - ~E[£(x + W t ) x W t ] dt, x G R 3 , 

2 Jo ^ 

where {W t )t>o is a standard 3D-Brownian motion. 

Proof. The probabilistic formula derives from Proposition 14.21 and the regularity of u from 
Propositions 14.21 and 14.31 We prove the uniqueness of the representation: since divw = 0, 
we have u = curl^, where ift is the stream function. Now, by the maximum principle, the 
unique solution of the problem 

Am = 0, u — > as |x| — > oo 

is u = 0. □ 

Since we are interested in the time evolution of the vector fields, it is appropriate to give a 
time-dependent version of the previous theorem. We recall that the spaces U a (T) and ZY^-(T) 
have been defined in <ED and (Q, the spaces V a ' p (T) and V£* P (T) have been defined in (Q 
and (EH). 

Corollary 4.5. Let a G (0, 1) and 1 < p < §. The map & : V a ' p (T) -> W Q (T), defined as 

Bs(£)tb,x) = \- r-E[at,x+W s )xW s ]ds, 
1 Jo s 

is linear bounded and ||£fc|| < C, where C is the constant, independent ofT, appearing in 
Proposition \4.3\ 

Moreover, ifL, M > are constant such that M > CL, then the map Gs : V^' P (T) — > U^(T) 
is linear bounded. 
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5. The representation map 
The section is devoted to the study of the properties of the representation map M, defined as 

M(u) (t, x) = E[U^(Xf)) + f E[U?g{t - s, )] ds, 

Jo 

where tp = ip(x), g = g(t, x) and Xp l are the Lagrangian paths, defined in (I5.3I) . and U x,t are 
the deformation matrices, defined in (I5.4I) . 

In the first part, some regularity properties of the Lagrangian paths and of the deformation 
matrices are obtained. In the second part we show that M maps the space U a (T) in V a ' p (T) 
(for the definition of the spaces, see (12.11) and (I2.2!) ). Finally, in the third part, we prove that M 
is Lipschitz-continuous from U a (T) to V a ' p (T). 

5.1. Regularity of the Lagrangian paths. In this section we study some regularity properties 
of the Lagrangian paths 

f dX* = u(s, X x s ) ds + V2vdW s , s G [0, T], 

I xs = x, 

and of the deformation matrices 

/ dUJ = U;V{s,XZ)ds, se[0,T], 

I us = 

where u G C([0, T]; C^R 3 , R 3 )) and V G C([0, T]; C b Q (R 3 , R 3x3 )) are given. Notice that both 
equations have unique strong solutions. Hence, fixed a 3D Brownian motion ((W s ) s >o, (J r s )s>o) 
on the probability space J 7 , P), for each x G R 3 there is a process (Xf, Uf) s > that solves 
the corresponding equations, and the solution is unique up to indistinguishability. The equations 
can be solved path- wise, choosing the u> G £1 for which s — > W s (uS) is a continuous function. 
Hence, the statements of this section are true for all such cus, independently of x and s. First 
define 

HU,« = su p ll i; ( r )lloo- 

0<r<s 

Lemma 5.1. Assume u G C([0, T]; C^R 3 , R 3 )). Then 

(5.1) \X X S -X y s \ < \x-y\e s ^ u ^% s > 0, x, y G R 3 . 
Moreover, if div u = 0, then for all s > and u G fi, map 

iGR 3 h G R 3 

is a diffeomorphism, the determinant of its Jacobian is everywhere equal to 1 and 

(5.2) [ (p(Xf(w))dx= f v{x)dx tp G L^R 3 ). 
Proof. First we prove (15.11) . By easy computations, 

\X X S - X y s \ <\ x -y\+ \\Vu\\^ s [ \X x r - X y \ dr 

Jo 

and applying Gronwall's lemma, we can conclude. 

Using Theorem 4.6.5 of Kunita lfT8ll (actually the assumption of Holder continuity on u is 
useless for our aim, since we deal with an additive noise, see also Theorem 4.1.1 of Busnello 
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0) one can easily deduce that x i— > Xf is a diffeomorphism and the determinant of its Jacobian 
is constant. Moreover, since div u = 0, the determinant of its Jacobian is equal to 1 for all times, 
so that, by a change of variables and a density argument, also (15.21) can be deduced. □ 

Lemma 5.2. Assume u G C([0,T]; C^R 3 ,!*, 3 )) andV G C([0,T]; C 6 Q (R 3 , R 3x3 )). Then 

\U X S \< e'lPH— x G R 3 , s G [0,T] 
and for x, y G R 3 and s G [0, T], 



a, ,s 



x - y 



« 



Proof. The first property derives from (13.81) . As regards the second, from (13.81) and (15.11) we get 

rS rS 

\U?-U?\ < / ||X>(r)|| 0O |E^-£^|dr+ / e sm ^ s \V{r,Xf) - V{r,X y r )\ dr 
Jo Jo 

< ||P|U, S / \U X -U^\dr + se s ^- +s ^ u ^[V] a , s \x-y\ a 
Jo 

and, by Gronwall's lemma, the claim follows. □ 

Let £>&(R 3 , R 3 ) be the space of all bounded Borel-measurable functions and define the oper- 
ator 

g^(x)=E[[/>(Xf)], x G R 3 . 

Lemma 5.3. Let s > 0, ?/zen 

1. Q s G £(£ fe (R 3 ,R 3 )) anJ ||Q s |U (Bb) < e'lPH— 

2. Q s G £(Q*(R 3 , R 3 )) am/ ||Q a |U(c«) < e 2 -" 27 "— + 'H V «H—(1 + a[P] a> .) 
Moreover, if div u = 0, 

3. G £(Lf(R 3 , R 3 )) and ||Q s |U (iP) < e'" 7 '"— 

Proof. First property is an obvious consequence of the previous lemma. About the second, 
using the two lemmas above, 

\E[u;<p(x;)-u*<p(xy)]\ < e\u: - uv\ ■ |^)| +E|t^| • - ^)l 

<-" fcTTIl J_ 1 ^a2s|IXJ||oo.»+s|IVw|[oo,a II ,-,11 I «. ,,\a 

Finally, assume div u = 0. Using (15.21) . Holder inequality and the previous lemma, we get 

/ \Q.<p(x)\ p < e HI^ E / MX?)| p <e^HMI£- 

□ 

5.2. Definition of the representation map. Here we prove that A£ maps U a {T) in V a,p (T). 

Before proving such claim, we need some preliminary definitions and results. For each u G 
U a (T), consider for all x G R 3 and t G [0, T] the Lagrangian paths 

dXf'* = -u(t - s, X'f) ds + V2vdW s , s G [0, t], 

yX,t 

Aq — X, 

and the deformation matrices 

dUf = UfV u (t - s, Xf) ds, s G [0, T] , 

JX, 

'o 



(5.3) 



(5.4) 
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where V u = Vu or V u = \{Vu + Vu T ). 

Lemma 5.4. Let u e U a {T) and if) e C 6 (R 3 , R 3 ) n L P (R 3 , R 3 ). The function 

(s,t)e{0<s<t<T}^ E[C/ S 'V(^'*)] e L P (K 3 , R 3 ) n C 6 (R 3 , R 3 ) 

is continuous with respect to both variables. 

Proof. First we show the continuity in C&. If < s < t < T and < r < v < T, with t < v, 
then for each x G R 3 , 

(5.5) < E|C/f - C/^l • + E|C/ri " llW) - 

In order to estimate the different terms of the above inequality, we see that from the equations 

\U? V - Ur\ = I jT U x /V u {y - a,XZ> v ) da\ < e^HlPJU^r - s\, 

and 

(5.6) \X? V - X*' v \ < \\u\U\s -r\ + V2^\W r - W s \. 
Moreover 

\X?-X? V \ < [ S \u(t-a,X x /)-u(v-a,X^)\da 
Jo 

f \\u(t -a)- u(v - ^lloo + f \\Vu(v - ^lUlXf - X? v \ 
Jo Jo 



< 

and, by Gronwall's lemma, 

I vx.t 

Finally, 

\U?-Uf' v \ < [ \U*/\-\V u {t-a,X*/)-V u {t-a,X^)\da 
Jo 



X*' v \ < e vM °°'" [ \\u(t-a) -u(v- a)^ da. 
Jo 



+ f W x /\ • \V u (t - a, X*f) - Vu (v - a, X^)\ da 
Jo 

+ f \V u (v - a, X? v )\ ■ \Uf - U x f\ da 
Jo 

< e^ll—H^lloo^suplXf -X*f\ 

[0,v] 

+ve v\Wu\u v gup ||x>„(« - a) - V u (v - ^Hoo 

[0,v] 

+ ||V«||oo,« \ S \V X / -V x f\da 
Jo 

< A(t,v) + C f \U X /-U^\da, 

Jo 
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where EA(t, v) — > as \t — v\ — > 0, and by Gronwall's lemma, 

\U X / - U*' v \ <A(t,v)e Cv . 

Using the above estimates in (15.51) . it is easy to show continuity with values in C&. In order to 
show continuity in L p , we remark that the above estimates ensure convergence for all x G R 3 , 
so that we need only to show uniform integrability. To this end, notice that, by a change of 
variables, 



/ \E[Uf^(XfW < CE [ my)\ p dy 

J\x\>K JX/(B C K ) 



< c f \^m\xf- x \<^) 

J\v\>f 2 

+C|MgP[|*?'* - X\ > y], 

< c I \^ + cu\\ p P nxf-A> 

J\v\>4 



where C = Te T ^ u ^°°- T , and, because of (15.61) . for K — > oo, the above quantity converges to 
independently of s, t. □ 

Now it is possible to prove the above mentioned result on the map M. 

Proposition 5.5. Given 1 < p < § and < a < 1, let ip G C fe Q (R 3 , R 3 ) n L P (R 3 , R 3 ) and 
g G V a ' p (T), then M maps U a {T) in V°' P (T) and 

(5.7) \\M{u){t)\\ c ^Lv <e 3 *ll v '"'l-'(l + t||V M ||c r )(||^l|c r nL P + / \\g(s)\\ c ^ LP ds) 

Jo 

Proof. First, M(u) G C£ D L p follows by Lemma 1531 moreover also estimate (15.71) can be 
easily deduced. Finally, from the previous lemma it follows that 

t^M(u)(t) g C£ n L p 

is continuous. □ 

5.3. Lipschitz continuity of the representation map. Letg G V a,p (T) and?/; G C^(R 3 ,R 3 )fi 
L P (R 3 , R 3 ), and consider the map 

M:U a (T) — ► V a ' p (T) 

defined in the previous section. The aim of the present section is to show that such map is locally 
Lipschitz-continuous. In order to do this, we will use Girsanov formula. First we rewrite M in 
a more appropriate form, namely 

M{u)(t,x) = E[F tjU (X x > t >*)], 

where for each trajectory w G C([0, T]; R 3 ), 

F t , u (w) = V{' u (w)ip(w t ) + / V s t ' u (w)g(t-s,w s )ds 
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and V t,u (w) is the solution of the following differential equation 



V s t > u = V s t > u V u (t-s,w s ), s<t, 
V*> u (w) = I. 



Notice that Uf' t,u (uj) = Vg' u (X x ' t ' u (u;)), for each lu e and we have made an explicit refer- 
ence to the dependence from u in the Lagrangian paths and in the deformation matrices 

JJX,t,U 

By Girsanov formula, we have 

E[F tiU pT^)] = E[Zf' u F t , u (x + VtoW.)], 

where 

= exp[—L f (u(t-r,x + V2^W r ),dW r ) - f \u(t - r, x + V2^W r )\ 2 dr] , 
V2u Jo 4^ J 

with s < t, so that for each «, 

A£(w)(t,x) = E[Z?' t > u V t t ' u (x + VtoW.)i/>(x + VtoW t )] + 

+ / E[Zf' u V*' u (x+ V2uW.)g(t- s,x + y/2uW t )]ds. 
Jo 

Using this representation, we will prove the following proposition. 

Proposition 5.6. Given 1 < p < § and < a < I, let if) G Q*(R 3 , R 3 ) n L P {K 3 , R 3 ) and 
g eV a ' p (T)andset 

£ o = \\i(>\\c?nLP + / \\g( s )\\c«nLP ds. 
Jo 

For each u, v G U^iT), 

sup \\M(u) -M(v)\\c°riLP < C(v,p)C M (T)e sup ||u - v\\ c i, a , 

t<T t<T b 

where C(u, p) is a constant depending only on p and v and lim^ Cm(T) = 0. 

The proof of the above proposition will be carried on using the subsequent lemmas. In order 
to make the explanations easier, we introduce the following notations. We define A xy f = 
f(x) — f(y) for any function /. Notice that 

(5.8) A xy (fg) = (A xy f)g(x) + f(y)(A xy g). 

If the functions depends on two variables, we define A uvxy as A uv A xy and, by applying twice 
the above formula, 

A uvxy (fg) = A uv [(A xy f)g(-,x) + f(-,y)(A xy g)} 

(5.9) = (A uvxy f)g(u,x) + [A xy f(v)}[A uv g(x)} 

+[A uv f(y)][A xy g(u)} + f(v,y)(A uvxy g). 
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Lemma 5.7. Let u, v G U^ f (T), then for each w, w' G C([0, T]; R 3 ) and for all s <t, 

IK*'" HI < e<M , 

\A UV V*''H\ <te 2tM \\u-v\y b , 

\A ww ,V s t ' u (-)\<2Mte 2tM \\w-w'\\ a O0 , 

|A„wK*''(-)l < (1 + 3tM)te 3tM \\w - w'\\^\\u - v\\ c i, a . 



Proof. The proofs of these properties are similar, we just give the proof of the last one. Indeed, 
using formula ( 15.91) . 

-^(A uW V s *'-(-)) = A uW (^^--(-)) = A nW (\A- (-)P. (t-s,-)) 

= [A uvww rV^(-)]V u (t - s,w„) + V*' v (w')[A ww ,V u _ v (t - s, ■)] 
+ [A ww ,V s t ' v (-)]V u _ v (t - s,w s ) + [A U X'~ 0')H A ww ,V u (t - s, •)] 
so that, by using the other inequalities of this lemma, 

|A uW ^-'(-)| < M \A uvwwl V^(-)\dr+\\w-w'\\^\\u-v\\ c i !a \V r ^(w')\dr 

Jo b Jo 

rS rS 

+||«-v|| c i, a / \A wwl V r ^(-)\dr + M\\w-w'\\^ / lA^^'-^Ol* 
6 Jo Jo 



< M [ \A uvww ,V r t '-(-)\dr+(l + 3tM)se 2 
Jo 

and, by the Gronwall's lemma, the inequality follows. □ 

Using the previous lemma and formulas (15.81) and (15.91) . we can easily deduce similar prop- 
erties for the functional F. 

Lemma 5.8. Let u,ve Um(T)> then for each w, w' G C([0, T]; R 3 ), and for all t G [0, T], 



F t ,uH\ <e' M [|^MI+ / \g(t-s,w s )\ds] 



\A uv F t ,(w)\ <te 2tM \\u-v\\ c iU(w t )\+ [ \g(t-s,w s )\ds] 

Jo 

\A ww ,F t>u (-)\ < (1 + 2tM)e 2tM e \\w - 

|A uW F t ,.(-)| < (2 + 3tM)te 3tM e \\w - w'\U\u - v\\ c i, a , 

where e = \\i/)\\c% + J * ll#( s )lk Q ds. 

Finally, we estimate the same quantities on the process Z. 
Lemma 5.9. Let u, v G U^iT) and q > 2. Then for all s <t, 

E|Zf 'T < Ce ctq/2M \ 
E\A uv Z^-\« < Ct^ 2 e CMqtq,2 \\u-v\\ q Cb , 
E\A xy Z^ u \ q < Ct q/2 M q e CMHq/2 \x - y\ aq , 
E\A uvxy Z^r < Ct 5q / 2 M 2q e CMqtq/2 \x-yn\u-v\\ q 



A PROBABILISTIC REPRESENTATION FOR THE VORTICITY... 31 

where C = C(q, v) is a constant depending only on q and v. 
Proof. From the definition, we see that Z* ,t,u solves 

dZf> u = -±=Z x /> u u(t - s, x + y/2uW a ) dW s , s < t, 

Again, the proofs of the four inequalities are similar, we prove only the last one. By applying 
formula (15.91) . we get 

ds{,A UVX yZ g ) A uvxy (d s Z s ) 

' \A uvxy Z/'-)u(t ~ s, 17) dW s + Zf' v A xy [u(t - s, Y-) - v(t - s, Y s )\ dW s 



V2z/ 

+{A xy Zf' v )[u{t - s, Y s x ) - v(t - s, Y s x )\ dW s + (A uv Z y /'-)[A xy u(t - s, Y s )\ dW s 



where, for the sake of briefness, we have set Y x = x + y2vW s . By the Burkholder, Davis and 
Gundy inequality, 



E| A uvxy Z/' \ q < C 



M q E[ / \A uvxy Z-/'\ 2 dr]i + \\u - v\\ q ca \x - y\ aq E[ / \Z y /' v \ 2 dr]i 
Jo b Jo 

+ \\u - v\\ q c E[ / \A xy Z^ v \ 2 dr]i + M q \x - y\ aq E[ / \A uv Zf-\ 2 dr]i 
b Jo Jo 

so that, by using the Holder inequality and the other inequalities of this lemma, we get 

E\A uvxy Z/'-\ q < CM q si" 1 [ E\A uvxy Z r ' t '-\ q dr + CM 2q t q s q/2 e CMH9/2 \x - 

Jo 



y\ aq \\u- v \\ q 

°6 



Finally, using the Gronwall's lemma, we obtain the required inequality. □ 

We are now able to prove the main result of this section. 

Proof of Proposition 15. 61 Let u, v E Uf^T). We start with the estimates in Cb and LP. Using 
formula (15.81) and Holder inequality, we get, for each x E R 3 and t > 0, 

\[A uv M(-)](t,x)\ = \E[A UV (Z^ F t ..{Y x ))}\ 
(5.10) < ^(^[(ElA^Zf'ry^lElF^^)! 9 ) 179 

+ (E\Zt ,t,v \ q ') 1/q ' (E\A uv F tr (Y x )\ q ) 1/q 

where q > 1, q' is the Holder conjugate exponent of q and we have set Y x = x + V2uW s . 
Using the estimates in Lemma l5~8l and in Lemma l5~9l and the inequality above with q = 2, we 
obtain the estimate in the Cb norm, 



sup||A u J\£(.)|k <Ce (T + VT)e( CM2+2M ^ ul , . 



U — V\\r,l,a . 



h 



t<T 

Using again Lemma l5~8l and l5~9l and the inequality (15.101) above, with q = we can obtain the 
estimate in the L p norm, 

sup \\A uv M(-)\\ p L p < Cel(Tv + T P/ 2 )e 2 ™ +CMP ' tP ' /2 \\u - v\\ c i.a. 

t<T b 
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To conclude the proof, we need the estimate in the C£ norm. For all x, y E R 3 and t > 0, by 
applying formula (15.91 ) we get 

\A uvxy M(-)(t,-)\ < E[A uvxy (Z^F t ,(Y))} 

< E[(A uvxy Z-/'-)F tiU (Y x ) + Zf ,t,v [A uvxy Ff r ( Y" )] 

+(A uv Zf'-)[A xy F t;U (Y-)} + (A xy Z t ' t ' v )[A uv F tr (Y x ) 
Using the inequalities in Lemma l5~8l and Lemma l5~9l it follows that 



\A uvxy M(-)(t, ■) | < Ce (Vt + T + MT Z ' 2 + MT 2 )e 3TM+c ™ 2 \x - y\ a \\u -v\\ ( , .. 



□ 
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